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Abstract 

Wc consider a class of random matching problems where the distance between two points 
has a probability law which, for a small distance I, goes like V . In the framework of the 
cavity method, in the limit of an infinite number of points, we derive equations for pk, the 
probability for some given point to be matched to its k-th nearest neighbor in the optimal 
configuration. These equations are solved in two limiting cases : r = — where we recover 
Pk = l/2 fc , as numerically conjectured by Houdayer et el. and recently rigorously proved 
by Aldous — and r — > +oo. For < r < +oo, we are not able to solve the equations 
analytically, but we compute the leading behavior of pk for large k. 

1 Introduction 

Let us consider a set of N points i = 1, . . . N, N even, and 'distances' between them Uj = A 
matching is a partition of this set into N/2 pairs. The 'length' of such a matching is 

L — ^ Ipair (1) 

pair&matching 

We focus on the properties of the minimal length matching. If the Uj are random variables, 
it is known that in many cases, the length of the optimal configuration and other quantities 
converge with probability 1 to their average value when N — * +oo. A lot of different probability 
distributions of the can be considered. In this note, we deal with the case where the lij are 
independent and identically distributed random variables, with a law p(l) defined on [0, +oo[. 
As noted in Q, in the infinite iV limit, the only relevant feature of p is its behavior around 0, 
say 



jr 

p{l) ~ - 



In the following we stick to the tradition that r is an integer, which simplifies results though 
it should not be difficult to generalize our computations to a non integer r. The thermodynamic 
limit of this model has been studied for a long time from the point of view of the statistical 
physics of disordered systems. The replica method yielded predictions for the mean length and 
the distribution of the lengths of occupied links in the optimal configuration Q . This was shown 
to be equivalent to a cavity approach ||. 
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Numerical works checked the validity of the results obtained with these techniques Q [|], 
and dealt with another quantity, the probability p^ for some given point to be connected to its 
k-th. nearest neighbor in the optimal matching. Houdayer et al. || conjectured that, in the case 
r = 0, one simply has 

Recently, Aldous || confirmed the mentioned predictions by rigorous proof in the case r = 0, 
and gave heuristic indications to generalize his method to an arbitrary r. 

In this paper, we show that insight on pf. can also be gained through a cavity approach, and 
in particular we recover (||). It gives further evidence that this method, as well as the replica 
one, exactly describes the matching problem, though on non rigorous grounds. 

The paper is organized as follows. In section ^, given any r we derive cavity equations for 
Pk and comment on their similarity with Aldous's ones ||. In section |3| we check that in the 
case r = we recover (||). In section |I| we derive, for < r < +oo, the leading behavior of pk 
for large k. In section |5| we compute pu in the limit r — ► +oo. Section ^ is the conclusion. All 
important results are supported by numerical simulations. 



2 The cavity prediction for 

The cavity equations at finite temperature have been derived in j|] (see also [0]). Let us briefly 
recall the essential steps — basically we reproduce ||] . One writes down a partition function for 
the matching problem 

Z= Yl ew(-f3N^L matching } (4) 

matchings 

l 

where the scaling factor iV r + 1 insures that the free energy F = — (1/(3) In Z is an extensive 
quantity. 

A tractable representation of this partition function can be borrowed from the theory of 
polymers Q. One introduces on each site a m dimensional vector spin Sj = (Sj, . . . S™) normal- 
ized by S| = m. Let d[i be the integration measure on the corresponding sphere. If we define 
T{j = exp ^ — /3-/V 1 / C^+i) ? one can check that 

Z = J™ / II Sjdfj,(Si) exp I ^ T^Sj.Sj J (5) 

m ^ » \i<3 ) 

thanks to the following property 

hm| d (i (8 i )S?S?>...S?=0 (6) 

unless p = 2 and ot\ = 02; in this case we have lim m ^o / ^/ i (S«)(5'f ) 2 = 1. The spontaneous 
magnetization of spin i has components mf = 6 a \mi. 

The cavity method consists in adding a new point i = to a system of iV points i= 1, . . . N, 
The partition function of the new system is computed assuming that the correlations of the TV" 
old spins can be neglected (technically speaking we suppose there is a single pure state): each 
spin lives in an effective field hi. So one writes 
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. N I N N \ 

Z N+ i = / II s l d ^i) ex P 5Z h i S i + T OiSo-Sj + h So J 
i=0 Vi=l i=l / 



(7) 



This is easily computed, thanks to @, and one gets by differentiation of InZ with respect 
to ho 

m = — (8) 

where we used the fact that mi, the magnetization of site i before the addition of the new point, 
is nothing but 

Note that in principle N should be kept even, but it has no incidence on our results in the 
infinite N limit. The occupation number of link — i reads 

noi = moT 0i mi (9) 

We are interested in the zero temperature limit (3 — > oo, for in this limit only the optimal 
configuration contributes. The probability distribution of the magnetizations in this limit has 
already been derived for another optimization problem, the travelling salesman problem JlCj] . 
We follow the same lines. It is useful to set 

mi = e 13 ^ (10) 

Let us then define 

Xi = N^TTloi-ifi i = l,...N (11) 
At zero temperature, equation (||) reduces to 



ipo = min Yi ( 12 ) 

t=l,. ..N 



and (pi) yields 



noi = 8i,i min (13) 

where i m %n is the indice attaining the minimum in dl^). 

In the thermodynamic limit, the ipi are uncorrelated and all have the same probability law 
P(<p) over the distribution of links. The x% are a l so independent random variables. Their 
common probability law is 



-oo r+oo 







n( x ) = / dl &pp{l)P{<p)8(N*il-ip-x) (14) 



1 f +(x > ( I \ 
= / dlp[^)P{l-x) 

jy— Jo \N~ J 

i r+oo rr 
~ —J dl —P(l - x) thanks to (2) 

Equation (|l2]) yields a self-consistency equation for P(<p) 

,+oc s N-l 



/ r+oo wv-x 

P(ip) = NU(ip) (j duU(u)j (15) 
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Let us define an auxiliary function 
We have M%) = G'(<p), and 



•roo ir+l 

dl - —P{1 - ip) 



(r + 1) 



(16) 



duTL(u) 



1 - / duU(u) 

-oo 



1 - jjGiv) 



so that dla) reads 



P{<p) = G>{<p)e- G ^ 



Plugging this into (16) we get an integral equation for G 



{t + i f Y e- G ^dt 



(17) 



(18) 



(19) 



We now proceed to the computation of the probability pk- Point to be connected to its 
k-th. nearest neighbor in the minimum matching means that two conditions simultaneously hold: 

(i) is connected to some given point io (N possible choices, e.g. io = 1) 

(ii) there are k — 1 points (C^~_} 1 ~ N k ^ 1 /(k — 1)! possible choices, e.g. i = 2, . . . k ) whose 
distances to is smaller than l 0io . 

So the probability of this event reads 



NN 



k-l 



Pk 



(k-l)\ 



N_ r+oo r+oo 

TT / dipiP(<pi) / dloipQoi) 

: ,i-0O JO 



.1=1 

N 

n ®(ioi-hi) 

=k+l 





A' 



n e ( z oi-w 



.1=2 



flQfa-iPi + N&iki-loij) 

.i=2 



where the Heaviside functions Q inforce conditions (i) and (ii). 
Let us define P as the following primitive of P 



du P(u) = 1 — e 



-G(<p) 



(20) 



We can then perform the integration with respect to the (fi, i = 2,...,N in (E 
integrals with respects to the Zj, i = 2, . . . N factorize, 



(21) 
The 



jyfc r+oo r+oo 

Pk = (fc _ i)i J dipiP(tpi) dl 01 p(l 01 ) 



h)i 



dlp{l)P(tpx+N^{l-l QX )) 



k-l 



(22) 



n N-k 



Ulo 



dlp(l)P(tp 1 + N~(l - I 



01, 
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After some changes of variables, it reads 



Pk = Ti TTT / dipP[ip) —p\ — 

OO J ]\[r+l XN^ 1 



(k-1) 



]\[r+l \N r + 1 J 



1 /"+ 00 du l + u\ ~. , 



n k-1 



N-k 



(23) 



Nr+l JO Nr + 1 \N r + 1 



Simplification of the first line above is straightforward using (^) . The second line needs some 
little work. It can be written 



+°° / I \ 1 f +OQ ( I + u\ ~ 

dvp[v + — + — / du pi — [P(tp + u) - 1] 

o V N~ J jV^+r \N~ J 



-i N-k 



(24) 



As p ( v + — 1 -t- ) = p{v) + —^-p'(v) + • • • + jp^TwP {r+l) {v) + ■■■, the first term in the 



bracket reads 



N r 



I + -A-K+00) - p(0)] + • • • + rt!,„ ^ r) (+oo) - P (r) (0)] + . . . 
_/V~ (r + l)!iV 

= 1 + + --- + FTim[°- 1 ] + --- 
_ r +1 

~ (r + 1)!JV + "' 

whereas the second term is simplified once again using (|2|). We get that expression (| 
to 



(25) 



tends 



ex P I ~ - I + / _, [ p (y + «)-!] 



(r + 1)! 7o 



r! 



(26) 



and finally we find for pi- 



I r+oo r+oo jr 

Pk = ~ TTT / dipP(lf) / dl — 

o Jo r\ 



{k-\)\ 



du - ^-P[}p - u) 



o r\ 



k-1 



(27) 



/ l r+1 r+°° (I + u) r ~ N 
x exp -T^W + I du —^{ p (v + u ) ~ 1] 



For some computations we may find easier to use a different formulation of the same result. 
Let us define the generating function of the p^ 



+oo 

H(e) = J2 VkE k ~ X 
k=i 



(28) 



Then H{e) is given by 
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r+°° r+°° r I l r+1 r+°° (1 + uY ~ 

H(s) = j ^ d<pP&) j o dl - exp i-^jy + J q du ^-[P( V + u) - 1] (29) 

+ e/ o '^j^-«)) 

Let us comment on our solution. It involves as a basic ingredient the function G determined 
by equation (IS), the other quantities being expressed in terms of G. This equation is actually 



well known, for it also appears in the solution with replicas (see H, eq. (22)), up to normalization 
conventions. The advantage is that, within the cavity approach, G is given a simple probabilistic 
interpretation in terms of the occupation numbers. This enables us to easily compute elaborate 
quantities such as pk, which is hardly tractable with replicas. 

Aldous notices that his method also yields equation (19) (see section 5.1 of pi), but we would 



like to stress that his approach and ours actually meet at an earlier stage: indeed our essential 
equation is not ( |l9| ) but rather the couple (|l2|)-(|l3|); the rest, including the introduction of G, 
is mere calculus. Roughly speaking, Aldous's method consists in linking the N — ► +oo limit 
of the matching problem to a matching on an infinite tree. Let us concentrate on case r = 0. 
Such matching is constructed introducing i.i.d. r.v. X — related to the root — and pQ)j>i - 
related to its children — , and the points (£i)i>i of a Poisson process of rate 1 (see section 3.3 of 
H). These quantities satisfy 

X = min(& - Xi) (30) 

i>l 

and the root of the infinite tree is matched to its child attaining the minimum. 

This is strongly reminiscent of (|i~2|)-(|i~3|). More precisely, if one relabels the Zoi so that they 
be sorted in increasing order, and consistently relabel the cpi, the probability law of (iV7oi)i>i 
tends to a Poisson process of rate 1. Thus ( |l^ ) becomes identical to (|30|), provided we identify 
Xi and <pi. 

Indeed, if we also take into account Aldous's heuristic considerations, similarities can be 
found even before (|l2|) : in section 5.3 of ||, Aldous gives indications for a possible extension of 
his method to the finite temperature problem. Using the same kind of r.v. as in ([}(]) he writes 
down the equation 

/+oo \ - 1 

X=\£e-*Xij (31) 

which can be identified with the cavity equation (||). Of course, the same way we derive (12) 
from (||) using trick (|To|), Aldous observes that (j3l] ) is consistent with (|30|) in the A — ► +oo limit. 

Such formal similarities are striking but we are unable to account for them. Whether some 
further study can lead to more a rigorous basis for the cavity method is not clear, because we 
define the , (fi through the somewhat exotic procedure of taking the m — > limit of a m-spin 
model. 

3 The case r = 

For r = 0, equation ( |T^ ) is easily solved. One finds 

G((p) = ln(l + e v ) (32) 
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r = 1 


r = 2 


r = 3 


(1) 


(2) 


(1) 


(2) 


(1) 


(2) 


Pi 


0.52891(15) 


0.52788 


0.53853(15) 


0.53865 


0.54308(15) 


0.54320 


P2 


0.22790(15) 


0.22816 


0.21706(15) 


0.21789 


0.21118(15) 


0.21192 


P3 


0.11215(11) 


0.11239 


0.10634(11) 


0.10702 


0.10318(11) 


0.10393 


PA 


0.05875(8) 


0.05882 


0.05719(8) 


0.05742 


0.05628(8) 


0.05648 


P5 


0.03185(6) 


0.03188 


0.03234(6) 


0.03243 


0.03247(6) 


0.03261 


P6 


0.01753(5) 


0.01766 


0.01881(5) 


0.01895 


0.01946(5) 


0.01960 


Pi 


0.00986(4) 


0.00992 


0.01126(4) 


0.01135 


0.01208(4) 


0.01212 


P8 


0.00556(3) 


0.00564 


0.00689(3) 


0.00692 


0.00763(3) 


0.00767 


P9 


0.00317(2) 


0.00323 


0.00424(2) 


0.00429 


0.00491(3) 


0.00494 


PlO 


0.00179(2) 


0.00186 


0.00264(2) 


0.00269 


0.00320(2) 


0.00323 



+oo 



(2) 
0.55523 



0.19329 



0.09399 



0.05278 



0.03222 



0.02079 



0.01396 



0.00967 



0.00687 



0.00497 



Table 1: pk computed (1) by averaging over 40000 samples of size N=200, (2) by numerical 
integration of our equations 



It turns out that the computation of H (e) is simpler than the direct computation of p^- Here 
|29|) reduces to 

/+oo / /-+oo \ >-+oo / rl \ 

dipP(p) exp (J du[P((p + u) - 1] J J die' 1 exp I ej duP(ip - u) J (33) 

r+oo e if ^ ,+oo / 1+e v 

= / dip— -K- / die exp I e In — 

7-oo {1 + e<*) 2 1 + e<t Jo P V 1 



+ e^- 



o 



+oo y rl fl + V_\ e 

V {l + yfJo X \l + yx 
+°° d y l + y-(l + yY 



(l + y)3 (l-e)y 
Provided e < 2 we finally get 



H(e) = — (34) 



which immediately yields the expected result (H). 



4 The case < r < +oo 



In this case, we are not able to find an explicit solution of (^). To check the validity of our 
approach, we solved this equation numerically by iterations for r = 1,2,3 and computed p^ using 
fl27|). The results compare very well with the estimates of p^ we got numerically by averaging 



over 40000 random instances of size N = 200 (see table |). 

Still we wanted to see if one could state as a simple law as (|3|), at least asymptotically for 
large k. Rewrite (f27|) as 

1 r+oo r+oo jr 

1 dp dlG'(v)-e A(k '^ l) (35) 



where 
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A(k,<p,l) = -G(<p) 



l r+1 f+°° , (l + u) 

VTTy. + L du — 



P{p + u) - 1 



(36) 



+fcln 



o 



1 du - 1 U)r p(<p-u) 



r 



To get the asymptotic behavior of such an integral, we focus on the /c-dependent maximum 
of A(k, ip, I) with respect to (tp, I). In the following we exhibit a stationary point and, assuming 
it is the absolute maximum, carry out the computation of the leading behavior of p^- 

We start deriving an asymptotic expansion of A(k, ip, I) for large k, p and I, under the 
constraints 



1< Ci < l/tp <C 2 < +oo (37) 
ip > k v (38) 

for some arbitrary but fixed C\,Ci and v > 0. We now deal in turn with each term in (|36|). 



If we write (16) as follows 



6Vi - / "dt ^f^ PU) 



(r + 1)! 

then expand (t + p) r+1 , and split the integral, we get 



(39) 



r+1 m r+oo j.r+1— m r+1 m p—tp j.r+1—m 

G{ip) = V r- / dt-—- -P(t) - y / dt-—- -P{t) (40) 

_ n m! 7-oo (r + l-m)! ml J-^ (r + l-m)! w v ' 



The terms of the second sum above all go to (it is easy to derive from ([II]) that P goes 
to zero faster than any power both in +oo and — oo). So 

( r + 1 ) ! ^0 

where the coefficients a m are essentially moments of the distribution P 

I r+oo +r+l—m 

a m = —. dt—- -P{t) (42) 

m! j-oo (r + 1 — ray. 



Using the definition ( |2l| ) of P, the third term in ( pq) can be written 

r f/ ,.\m /-+oo 4.r-m 

m! ^ (r-m)! 

All terms go to zero — it is essential here to be able to compare p and I; that is why we 
set®. 



8 



Finally, the argument of the In in the fourth term of (|3^) can be integrated by parts, giving 



k In 
k In 
k In 



r +1 

(r + l)! 

l r+1 
(r + l)! 



1 - e- G ^)" 









-i (r + l)! 

+ - { — L -Eg- ^ m °- + °(!) 



(r + l)! 



m=0 



(44) 
(45) 



(r + l)! r + l)! 



m=0 



+ o(l) 



because the o(l) terms at line (45) go to faster than any power of k thanks to (|37|)-(|38|). 
Collecting all the pieces, we get 



A(k,<p,l) 



+kln 



jr+l 



+ (' + 1)1 h 



r+l 



a m (l - ip) r 

m=0 



(r + l)! 

Dropping the o(l) term, the stationarity equations read 



+ o(l) 



(46) 



^(m+l)a m+ i(^ m + A:- 

m=0 



r-1 



+ £L=o( m + j>m+l(j ~ l) T 

E^ =0 « m (/-^ m 



(r+l)! 



r-a-yV; ^r-i 



+ £:- 



r\ 



(r+l)! 



Em=0 Om(! ~ ¥>)' 



(47) 



Let us call (fk,lk) the solution of this system. We can see that both ip k and l k scale as 
£,i/(r+i)^ w j 1 j c ] :i i s consistent with hypotheses (|37|) and (|38|). If we write 



h 



k^L 



the system reads 



(48) 



r-l 



-— — > k r + l 



T 



m=0 



+- 



(r+l)! 



-Em=o fcr+1 a m (l k -ip k ) r 
K , - EL=o k^(m + l) Qm+1 (Z fc - y fc ) m 



I + 



(r+l)! 



J2 r m =o kr+1 a m (i k -(f k y 



(49) 
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We face a perturbative problem whose small parameter is k 1 /( r+1 ). Knowing the a m , it is 
in principle possible to compute iteratively the following asymptotic expansions of p>k and 4 

r+l 

Pk = fioo + 2^ c m k -+ 1 + o(k~ ) (50) 

m=l 
r+l 

4 = ko + 2_, d mk r+1 +0(k~ 1 ) 
m=l 

Though we cannot compute all the coefficients in closed form for a generic r, we can give a 
simple procedure to determine the leading ones, loo and poo- We define a = loo/ Poo- Combining 
the two equations (^), where one sets k = +oo, yields 

(a - l) r ~~o?~ x ^ 

Note that a must be greater than 1 to be consistent with hypothesis ([3?]). On ]l,+oo[ the 



left hand side of (51) is a decreasing function of a, from +oo down to 0, so the above equation 
does admit a solution in this interval. Other interesting properties of a are derived in appendix 
A. The first line of system ( |49| ) gives 

« 1 = ( r + 1 ' ! ^TT (52> 

from which follow two relations we will need further down 



~r+l , tr+1 



V 



r+l 



(J 



oo Poo 



\r+l 



(r + l)! 

( r + l)!(a-iy 



(53) 



In ( |35j ) we perform the rescaling 99 <— p^p and i <— so that 



1 f 



r+l 



Pk+1 ~ ^|^^T 



+00 



dp / dll r G'(p k p)e A ( k '^M 
Jo 



(54) 



To get the asymptotic behavior of the integral, we expand A(k, p&p, l k l) around its maximum 
(p = 1, 1 = 1), up to quadratic terms, and extend the range of integration with respect to I down 
to —00 



1 c 1 



2tt 



«! H VdetdM 



G"(p fe )e A(fe '^ A) 



(55) 



where detd 2 A = ^yO^yi) _ ( ^WAO ) 2 |^ =1 . Looking at (p and @ it is 

clear that (|50| ) induces an expansion of A(k, pkP, hi) in powers of /c 1 /( r+1 ) of the following form 
(for finite p and I now) 



A(A W ,+. /,/./,•) -k Mr+ * + Jf ~ v)r+ \ £ k^fm(pJ) 

r + l)! ^— ' 

m=0 



(56) 



+Hn 



(r + l)! 



m=0 



+ o(l) 
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In this, ln[- • •] can be further expanded, yielding 



A(k, l k l) = _ k ^ l ) r+1 + ^) r+l +klnk (57) 

(r + 1)! 



+kln 



(r + 1)! 



r+1 



+ ^ k r + 1 Sm (<p, l) +0(1) 

m=0 



It follows that det<9 2 ^4 oc k 2 , (the symbol oc stands for 'goes like ... up to some constant 
factor'). The prefactor is a complicated but easily computed function of a, which we omit to 
write. 

Equation ( |4l| ) tells us that G'(<pk) ~ <Pk/ r ^ so that the factor <Pklk +1 G'((fk) in ( |55[) is oc k 2 . 
Beside, using (|53|) and (p7l) we get 



A(k,ipk,k) = -k + k\n.k + kr1n(a- 1) + ^ X m k^ + o(l) (58) 

m=0 

where A m = s m (l, 1). 

Dropping numerical prefactors, ( |55| ) thus reads 



^ke- k k k [(a - Iff exp ( £ A m ^) 

\m=l / 



Pfe+i oc -fce- fc fc'' [(o - 1ST exp ( X A,,,/." ) (59) 
Sterling formula finally yields 

Pfc+ icxVfc[(a-l) r ] fc exp( Am^j (60) 

\m=l / 

We have above the leading behavior of Pk+i- It is clear that an identical formula holds for 
Pk, but with different constants X m . 

This is certainly less appealing than the simplicity of case r = 0. The leading behavior of 
Pk involves more and more elementary functions as r increases. It makes us understand why 
Houdayer et aJ. || were not able to find any convincing fit of their numerical estimates of pk 
for r > 0. 

Note that as a consequence of ([30|), we can state the simpler but weaker following property 

lE£* _> r w a _i) (61) 
k 

This limit is negative, as expected, increases with r and goes to like l/2 r when r —* +oo 
(see appendix A for details). 

To check prediction ( |60|) we carried out numerically the whole computation sketched above, 
including constants, in the particular case r = 1. We got 

p k ~ 0.3387Vfc(0.6180) /c e- a6014 ^ (62) 
It is compatible with our numerical simulations (see fig. Q). 
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Figure 1: In the case r = 1, plot of the ratio [asymptotic behavior of equation (|62})1 / [numerical 
estimate of pt by averaging over 150000 samples of size N = 300] versus k 

5 The case r — > +00 

Our aim is to compute pk in the limit r — * +00, where new simplifications appear. We get 
inspiration from a similar computation performed in ||. We define the following function 



G(x) = G 



1 x\ . ,._JL_ 
- + ~ (rr\)r+i 

2 r J 



In terms of G, (1l9|) reads 



G(x) 



+00 / 

1+ - 

-r— x \ 1 



x + y \ r 



-G(y) 



dy 



When we send r — > +00 this equation has the following limit 

r+00 



G(x) 



*+y e -G(y) dy 



(63) 



(64) 



(65) 



whose solution is very simple 



(66) 



G{x) = e x 

As in the case r = 0, our choice is to first compute H{e). It is easy to explicitly rewrite (|2 
in terms of G, using (O) and (ET1) 



H(s) 



+00 



(67) 



x exp 



(e-1) ( {U + ^T- Tdt^Le-^ 



(r + 1)! 
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Then the changes of variables t = (1/2 + v /r){rr\) l ^ r+l \ u = (1/2 + w /r)(rr\) l ^ r+l ^ and 
ip = (1/2 + tp/r){rr\) 1 l { - r+ ^ give 



H{e) 



+00 - f + OO 

dtpG'(^)e- GW / dw e^+Ve-GH 



x exp 



(e - 1) e™^ - / dve v+w e- G ^ 



+oo 



dip- 



(1 - e) + - 1 



+ 1 



h (1 



e) [t + e-* - 1] + 1 
Remembering (pq), this means that 

-°° t dt 



1 



fc-i 



We are not able to do any further analytical simplification. However we can derive the 
asymptotic behavior of The integral in (p9|) is dominated by the region i — ► +oo, so that 



hoo / 1 \ fc— 1 



/>/, ~ / | tfte-Ml-^J (70) 

+oo _ / 1 \ fc— 1 



v 7 ^ / °° due-^ u 



The argument of the second exponential above has its maximum in u = 1. A saddle-point 
method yields 



y|fcV4 e - 2 ^ (71) 



The values of p& obtained by numerical integration of (|69D for r = 1 , . . . 10 are given in table 
|l| Extrapolating upon comparison with results for r = 0, 1,2,3, one can see that for small k, 
say up to k ~ 6, pj~ little varies when r increases. And these are the values of k for which p/, is 
significant. For larger values of k the range of relative variation is more important, and diverges 
when k — * +oo (compare (|7l| ) and (||)). But then anyway p^ remains very small. Beside, it 
seems that p2, P3 and p^ decrease with r whereas all other pk increase. 

We have not tried to see whether some of these features hold for the Euclidean random 
matching problem in dimension d — where the lij are the usual distances between N points 
independently and uniformly distributed in the hypercube [0, l] d , which is believed to be equiva- 
lent to the uncorrelated version we study here in the limit r = d— 1 —* +oo []|, 11]. In particular, 
one expects pk in this limit to be the same for both versions. 
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6 Conclusion 



We have shown that the cavity method yields predictions for which are consistent with both 
rigorous results and numerical simulations. 

Let us emphasize that our results also hold for the bipartite minimum matching (or assign- 
ment) problem. It has already been observed that the latter and the problem we discussed in 
this paper are equivalent in the N — > +00 limit, up to normalization conventions [||, In 
appendix B we briefly show how this can be seen within our formalism. 



Appendix A: Derivation of some properties of a 

Here we write a r for a, to highlight its dependence on r. Let f r be the function given on ]1, +00 [ 
by 

f r {x) = 1 - -1 (72) 



f T strictly decreases from +00 down to 0, and a r is defined by 

f r (a r ) = 1 (73) 

A small computation leads to 

M2-^) = ! + ^ + °(^) ( 74 ) 



and 



so that, for large enough r, 



Beside, 



> 1 for large enough r 



/r(2) = 1 - Yr < 1 ( ?5 ) 



2 " < ar < 2 (76) 



/r . +lK Ti) = ^_i J_ (77) 

( a ;+i _ l)r+l «r 

r r 

One has a r r+1 < a r and so (a r r+1 — l) r+1 < (a r — l) r+1 < (a r — l) r thanks to the right hand 
side of (f76|). Thus 

/ r+ i(a^) > / r (a r ) = 1 (78) 

It follows that 
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< "r+l 



(79) 



i.e. 



1 + — > 1 + 



(a r 



a' 



1) 



n 



r+l 
r+l 



> 



(q 



r+l 



by (73) 



(a r — l) r increases with r, and goes to 1 when r — > +oo, as shown by (|76|). More precisely, 
it goes exponentially fast to this limit: (|7^) implies 



1 



i.e. 



(ar - l) r 

(ar " 1) 



a, 

1 



+ o(— ) 
2 r 2 



by (76) 



(80) 



Appendix B: Note on the bipartite minimum matching problem 

The bipartite minimum matching problem can be stated as follows: take 2N points i = 1, . . . 2N 
divided into two subsets, A (i = 1, . . . N) and B (i = N + 1, . . . 2N) such that the distances 
between a point in A and a point in B follow the probability distribution p of (^), and allow 
only matchings where each pair is made of a point from A and a point from B. One can see 
it as a particular case of non bipartite matching problem for the points in iUB, setting that 
the distances inside each subset are infinite, i.e. Tij = for £ {l,...iV} 2 or € 

{N + 1, . . . 2N} 2 . Then imagine that you add a new point i = to, say, subset A .The cavity 
equation (S) is slightly modified 



1 



m 



J2i=N+l ToiTTli 



(81) 



where 



-l3(2N)r+H 0i 



(82) 



Provided we relabel the rrii, i = N+l, . . . 2N and redefine f3, we get exactly the same equation 
as in the case of the non bipartite matching problem. Indeed it is known from the replica 



approach [12] that, for systems where the lengths lij are uncorrelated random variables, the 
bipartite and non bipartite matching problems have the same limit when N — > +oo and differ 
in the 1/N corrections. 



References 

[1] J. Vannimenus, M. Mezard. On the statistical mechanics of optimization problems of the 
travelling salesman problem type. J. Physique Lett., 45 (1984) L-1145. 

[2] M. Mezard, G. Parisi. Replicas and optimization. J. Physique. Lett., 46 (1985) L-771. 

[3] M. Mezard, G. Parisi. Mean field equations for the matching and the travelling salesman 
problems. Europhys. Lett, 2, 913-918, 1986. 



15 



[4] R. Brunetti, W. Krauth, M. Mezard and G. Parisi. Extensive numerical simulations of 
weighted matchings: total length and distribution of links in the optimal solution. Euro- 
phys. Lett, 14 (4), 295-301, 1991. 

[5] J. Houdayer, J. H. Boutet de Monvel, and O. C. Martin. Comparing mean field and Euclidian 
matching problems. Eur. Phys. J. B 6, 383-393, 1998. 

[6] D. J. Aldous. The £(2) limit in the random assignment problem. math.PR/00W063. 

[7] G. Parisi. Spin glasses and optimization problems without replicas. Les Houches, Session 
XLVI, 1986 — Chance and Matter, ed. J. Souletie et al. (Elsevier Science Publishers B.V., 
1987), p. 525. 

[8] M. Mezard. Spin glasses and optimization. Invited paper at the Heidelberg colloquium on 
glassy dynamics and optimization, Lecture Notes in Physics 275 (Springer Verlag 1987), p. 
354. 

[9] P.-G. de Gennes. Exponents for the excluded volume problem as derived by the Wilson 
method. Phys. Lett, A38 (1972) 339. 

[10] W. Krauth, M. Mezard. The cavity method and the travelling salesman problem. Europhys. 
Lett, 8, 213-218, 1989. 

[11] M. Mezard, G. Parisi. The Euclidean matching problem. J. Phys. France 49 (1988) 2019- 
2025. 

[12] M. Mezard, G. Parisi. On the solution of the random link matching problems. J. Physique, 
48 (1987) 1451-1459. 



16 



